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PREFACE 


This investigation was prepared during the period 
February =- June 1952 at the United States Naval Postgraduate 
School, Monterey, California. 

The subject was suggested to the author by previous 
experience in maintenance of small high-speed diesel 
engines, during wnich burnt pistons had to be replaced 
on several occasions. 

The author wishes to extend his appreciation to 
Assistant Professor E. E. Drucker for his interest and 


advice throughout the development of the work. 
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The objective of this work was to establish by 
annalytical means the temperature distribution in an 
internal combustion engine piston. For this purpose the 
piston was divided into a disk and a barrel, for which 
boundary conditions were expressed for the exterior 
surfaces and a common conical surface. The problem was solved 
in two parts: first, a solution was found for the case where 
the gas temperature was assumed representable by an average 
temperature; this solution was then modified to account 
for a sinusoidal fluctuation of gas temperature. 

The resulting temperature distribution, within the 
limits of the approximations made, is: 
for the disk: 
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CRAPTER 1 
INTRODUCTION 


Upon analyzing a piston of an internal combustion 
engine with a viewpoint to understanding its behaviour 
as a heat dissipating member of the engine, it is readily 
realized that the piston absorbs a certain amount of heat 
through the upper surface of its crown. This heat is partly 
stored in the body of the piston until it reaches its steady 
state temperature; at the same time there is heat dissipation 
from the outside surface of the barrel and rings and also 
through its inside surfaces. After the steady state is 
reached, all heat absorbed mst be dissipated through these 
surfaces. 

According to 0. Le. Adams (1), five percent of the heat 
liberated by the combustion of the fuel mst pass through 
the piston crown surface into the piston. Of this heat, ten 
percent is dissipated by the lower surface of the crown of 
the piston, as stated by J. L. Hepworth (5). 

Radiation from the gases into the piston only takes 
place through a very small portion of the time cycle and 
it amounts, according to B. Pinkel (9), to only ten percent 
of the total heat being absorbed by all the metalic parts. 
forming the container for the hot gases; thus it may be 
seen that the heat absorbed by the piston itself through 
radiation will be a part of this ten percent depending on 
the relative areas of the piston crown with respect to the 
area of the rest of the combustion chamber. Throughout this 


work it will be assumed that all heat absorption by the 
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piston takes place by convection. Also all heat dissipated 
by the piston will be assumed to do so by a process of 
convection. 

The temperature of the gases throughout the cycle 
varies widely; the temperature curves for an engine may 
be obtained either experimentally by direct measurement or 
annalytically from knowledge of the cycle under which the 
engine is operating. In this work the temperature of the 
gases will be assumed, for the sake of simplicity, to 
follow a sinusoidal variation about an average temperature. 

The main objective of this work will be to determine 
the temperature distribution throughout the piston. For 
this purpose a set of boundary conditions will be prescribed 
wnich resemble as closely as possible the actual conditions 


under which the piston operates. 
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CHEPTER Il 
STATEMENT OF THE PROBLEM 


1. Description of piston. 

A piston with the characteristics shown in figure 1 
was selected for analysis. The crown will be a flat disk 
of uniform thickness 1. The barrel will be a cylindrical 
sleeve of constant thickness b. 

The boundary surface of disk and barrel will be 
surface a-a as shown in figure 1. This surface will have 
properties common to both disk and barrel and this will 
enable some of the unknown quantities entering the questions 
to be found. 

<~e Temperature scales. 

Figure 2 shows a sketch of the temperature pattern 
that a particular point of the disk is expected to follow. 
This figure is given mainly to show graphically some of the 
values being used throughout the development. 

Figure 3 serves the same purpose as figure 2 with 
relation to the barrel. 


The temperature of the gases will be assumed to be: 


Th ox = lle ~ WV. Am ATKT 


d 


where n will be the frequency of the cycle. 

The temperature Tg surrounding the outside surface of 
the barrel will be assumed constant. 

The temperatures for the disk will be measured from 


the level designated as in figure 2 and will be expressed 
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for @ny poimt of the disk. aims 
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The initial temperature will necessarily be a constant 
and will be designated as Bo. 

The temperatures for points of the barrel will be 
measured in the p scale of temperatures as shown in figure 


5, where: 


> = ¢ ae 


The initial temperature ¢, will be a constant. 
It is realized that the value of Tg will initially be 
equal to t,, but it is impossible with any degree of ease 
to consider into the problem its initial rise from t, to 
Tag and thus it will be assumed constant and equal to Tg. 
After a solution is reached, the results will be 
reverted to the standard level of measurement of temperatures 
in degrees F. which is designated as the t scale of 


temperatures. 


5.» Boundary conditions. 

The following set of boundary conditions has been 
selected: 

(a). At the lower surface of the disk, the rate of 
heat rejection is ten percent of the rate of heat absorption 
at the upper surface. Since the rate of heat conduction is 
a function of the temperature gradient normal to the surface 
considered, and if we assume the coefficient of heat conduction 


to be constant throughout the piston, we may state: 
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There will be a difference in the areas of the upper 
and lower surfaces determined by the relative values of 
Ry and Ry. In this work we will neglect this difference, 
but it is recognized that it could be included by mlti- 
plying the value given as 10 by the ratio of Ry to Ro 
squared. 

(b). At the upper surface of the disk the rate of 
heat absorption from the gases will equal the rate of heat 


conduction from the surface into the disk, or: 
L (2) = hs (50s ~ Ges) 
N4 ue 0 


(c). At the surface a-a the temperature of both the 
disk and the barrel will be equal. Since these temperatures 
will be measured from different levels, this relation may 
not be used directly without introducing undesirable 
constants. 

Therefore it is necessary to realize that at this 
boundary the isotherms from disk and barrel must coincide 
and furthermore they mst be continuous functions of the 
variables involved. This may be stated symbolically: 


along surface a-a;: 
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(d). The lower surface of the barrel will be assumed 


insulated, or: 


(Wh... *2 
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(e). The inner surface of the barrel will be assumed 
insulated, or: 


ean 


Jv ABR 


(f). At the outer surface of the barrel the rate of 
heat rejection from the barrel into the surrounding medium 
will be equal to the rate of heat conduction from the 


barrel into its outer surface, or: 


«(Bho (ae A) 


(zg). At time = zero, a point: 270, rz=0O, will be 
at a temperature o, or: 


at 
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CHAPTER IIL 
SOLUTION OF PROBLEM 


1. General solution. 

In solving the problem the disk and barrel will be 
treated separately and the solutions matched at the surface 
A-A. 

Both the disk and the barrel involve heat transfer in 
a cylinder; the disk may be considered as a flat solid 
cylinder, while the barrel may be considered as a hollow 
cylinder. Under these conditions the best solution of 
Fourier's Law of Conduction of Heat will be that expressed 


in eylindrical coordinates; 


die CCH commen 15 (1) 
Hea (HS A 1. 


Assuming the variables to be separable a solution as 
found in Appendix I iss 


ee (2) 


t 2£ e [Coby + Dain AT ereBh 0d] 
From figure 2 it may be seen that the expected curves 


for temperatures will be of the form: 
-aw 
o-oo, 2 ° 
for a particular point on the disk. 
At the same time a sinusoidal variation will be applied, 
which proceeding along the z-axis will be damped and out of 


phase; therefore our solution would have a form such as: 
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where cz is the phase angle and edz 


accounts for damping. 

Examining the previous equation it may be seen that 
the variables z and T may not be separated as assumed in 
the solution found by equation (2), and thus this equation 
will not yield an exact solution. 

A direct solution of the problem was attempted by the 
use of Duhamel's Theorem (3), but it was found impossible 
to match the result obtained for the disk with that for 
the barrel at the surface ae-a. 

Under these conditions it was decided to divide the 
problem into one where the gas temperature was assumed to 
be the average temperature Ty, and then proceed to supere- 
impose on it the effect of the sinusoidal variation in gas 


temperature. 


| 2 Solution for the disk. 


From equation (2): 
OE = +b "Te amb . D ain by // A So(ar) oO 


Since 9 mst be finite throughout the body: 
B= 0 
and therefore; 
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From equation (3): 


Jo ~a (a *4 IT ; i iz 
Jy > Fe I. en tun b 3 ¢ Dice b, ] SoCs) 
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From boundary condition (b): 
Jd _ = 
t ($F) a =i (6... O..) 
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The solution of this equation for b may only be done 
graphically; it will yield a series of values for b which 
will be designated as b,. 


Equation (4) now becomes: 


Oo 


— aay )T 
& : ds éa™ "fem bn 3 + pe anh EG (5) 


mri 


From boundary condition (g): 
T:0 (eee fu +o QO. F- 
Substituting this in equation (5): 


G: 8, 
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4- Solution for the barrel. 
From Equation (2): 


b . Geen tLeingy [[Rd (fu) +m K (fo) | (7) 


From boundary condition (d): 


{ 3) 0 
§ “get 
therefore: 


lia = tang bo 
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From boundary condition (e): 





at). 4 

V2 Jae eR, 

or: IY, Fj ({R) fF 
." 1 


therefore: 
Za ie tr 74) 
g- ae “h , ew 5 (te [BUI + F by) (9) 


From boundary condition (f): 


k (iH) == 4, (6,72) 
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Solving this equation graphically for the value of f, 


a set of values f, are found that satisfy this equation, 


p 
and equation (9) becomes: 


p : a oie (3-49) [eipu)% Ya) (10) 
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4. Conditions at the surface a-a. 


Recapitulating: 
@ >a = 
— aeaieee. ) —_ 
G : ?_ 90 e@ [em aa aint [Sead (6) 
WA 3] kh 
and : 
a: a am 
b = > eh +5 °)7 cans (s-Lo) Self. 2) + r, Kipadl (10) 
pei 
From boundary condition (c) at points along the 
surface a-a: J7O Jp 
J3* : Ja* (11) 
sso. J"? 
A Ae Jba™ (12) 
For equation (11) it is easily seen that for n=l 
and for a particular point along the surface a-as 
90 ces 
mT. =>) SS (13) 
“=| pi! 


then for n: 5:3 


oo Sf 
—= a eae 
2 Pn Lx. = Tolle: Sp tied 


pam 


this suggests that as a method of approach to the 


evaluation of a, g, and P we assume; 
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Acceptance of this assumption involves the necessity, 
in order to maintain the equality of the above equation, 
that: 

Ty, = Sy 
where: m=p 

Also, in order that the effect of time on the tempera- 
ture be equal for the same point on the surface a-a as 
approached by either equation (6) or (10), that: 

- & 


fy 


Applying this condition to equation (13), we find: 


6, [a A bag | # Pat bly nti api- SE fea ~ 


it he Vp RIT {pa) 


where: 
Hi = p 


This equation is valid for points along the surface | 
a-a2, were: ‘: < (&, smi) 
and should yield the values of Pp. 

Now, the replacement of equation (15) into equation 
(14) will yield an equation in terms of trigonometric and 
Bessel functions of r, valid for values of r ranging from 
Ri to Roe Since this equation may not be simplified to any 
appreciable extent, it has been considered more advisable 
to replace the values of z and r of a particular point 


along the surface a-a, such that: 


A/2 = No 


(Ry + Ro)/2 = Ra 


N 
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then we finds: 
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kK S.. Y. (YR) Tol, Re) (16) 


where: m= p 
If now in equation (12) for n = 1 this same point in 
the surface ae-a is investigated: 


pmb. - 4, orb b= pork. (t-Lilf- Ee Lf ta (17) 
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Equations (16) and (17) should give the same value 
for Poe There is no way of proving that this will be so 
except by working out a particular problem. If in so doing 
it is found that the values of Pp satisfy simltaneously 
both equations, this will then mean that the assumption made 
that: & = bp 
and the derived expressions: 

Tm = Sp 

and: fp = a 
are justified. 

In case that the values of Pp found by means of one 
of the equations mentioned do not satisfy the other 
equation, there is still the possibility of introducing 
this value in the second equation and determining a 
relationship between Ry, b, Lo, and 1l such that both 
equations are satisfied and for which the assumption made 
is valid. This would introduce a limitation to the physical 
proportions of the pistons for which this development is 


useful. The relationship just mentioned would take the form: 
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5. Partial solution using Tp. 


On the light of previous conclusions, the solution is 


now s 


ge ce —~ *, bi )r j 
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2 _2 weuiae, b )T ‘ 
d yD P, e€ io Cre b,.. (5-6 IIa) + Fe wh | (20) 
ee | ft! 
and since: i Tyas 
and: p ->C~ Ta 
we find, for the disk: 
2 ee (f+ 6, )& : 


and for the barrel: 
t: Ta + 


et / 


24 i ; 
2 _ < wh, pe (5-te)]. o( fpalt F e(e)] (22) 
‘aad 


6. Superposition of sinusoidal variation. 
In order to superimpose the effect of the sinusoidal 
fluctuation of the gas temperature about a mean temperature 


Tn» the above equation (21) in a simplified form: 


C = Toes Te ra CO a 
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after @ qwerter of a Cycle, or at T we find that: 


T gas =i +To Ren 2, & 
therefore: 
Tse = Ty, a Ole 


If this were put into equation (21) it would give an 
answer corresponding to the situation where this new 
temperature would have been imposed on the piston since 
time zero. 

To avoid this, instead of: 

Twi. Via T (23) 
use the following expression given by M. Jakob and G. A. 
Hawkins (7): | 
To e Abn (arnt = 31 ) 


(24 ) 


mere: 

T, is the amplitude of the oscillation at the surface 
of the disk; 
le accounts for damping effects along z, and 
- SE accounts for the change in phase along z. 

The above expression was given for a thick plate sub- 
jected to conditions similar to those of this problem. Even 
if in this case the disk may perhaps not be considered a 
thick plate, it is felt that this equation will yield a 
fairly good approximation if it is considered that it is 
common knowledge that the effects of fluctuating surface 
temperatures do not penetrate to any great extent beyond 
the surface exposed. In this particular case this effect 
is further minimized by the fact that the fluctuation of 


temperatures is very rapid. 
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From equations (23) and (24) it may be seen that what 
happens at the surface at time zero will have its effect 


at a depth z with a time lag of: 


} 
L\ (2 VE A Vain 


The addition of the term expressed in equation (24) 
will affect in equal form both the equations (21) and (22). 
Since for the barrel the values of z will, with the excep- 
tion of its topmost part, be far greater than those for 
the disk, this effect will only be perceptible in its 
upper section. 

Including equation (24) into equations (21) and (22) 
we have as a final result: 


; et 
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For ‘the disks: 


For the barrel: 
Th 


aa/= == 
ae | ae (avnT~27e ) + 


os? eth abe a - (L/S pad KU, | 


ees ps 

7. Conclusions. 

It 1s realized that the results obtained are not rigorous 
but it is felt that they constitute a fairly good approxi- 
mation to the actual temperature distribution in the piston. 

Further work might be done on the subject by solving 
an actual problem in order to check the practicability of 
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the method used to determine Pp. By actual experiment on 
an engine the overall results of this work could be tested 
for accuracy and the percentage of error, if any could be 


found. 
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APPENDIX I 
GENERAL SOLUTION 


To find a general solution to Fourier's Law of Heat 
Conduction for the case where the temperature should be 


analyzed in the unsteady state for a body best described 


in terms of cylindrical coordinates, we use; 


Eee, tO, UE 
IY = (ff+45, * jay (2) 
where; e 
> $( aa) T) 
Assume: { = O(+)* Rie) « ZG) 
then: 
— see 
hb AT KR dar AR dv 2 Ayr (2") 
Since r, z, and T are independent variables, we may 
say: 
ak Ak + ae AK 2 - a (3) 
Ant ne an 
— a - 
ee 
| dG 2 o> (5) 
a pea >c-_G_- 
dG Avr 


Equation (3) may be rearranged as: 


Mo. 2k: 0 
dr* NR dn 


which is Besselts equation of the first kind of order zero. 
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Then: R =AT(o1) +B Ke Car) (6) 


Equation (4) may be rearranged as: 


ge "7 = & 
A 3? 


which yields: 2> Cumb;+ D am by (7) 


Equation (5) may be rearranged as: 


3 + LQ (at+b ) = 0 
pa (at+b*)r 
which yields; @2€e (8) 


But, since: t = O~KRxzZ 


we find using equations (6), (7), and (8), that: 
wela*sb° IT /- 
s° Ee i ill + Dainb; [fF Toler) +B Yo(ar){ 
In this equation the constants a, b, A, B, C, D, and 
E mst be found so as to satisfy the boundary conditions 


of the specific problem involved. 
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